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Abstract 

The 3A^ forces due to tt — p,tt — a and tt — lo exchanges following from the nucleon Born diagrams 
and diagrams with an intermediate A^*(1440) are re-analyzed. The cancellation between vr — o" and 
vr — a; forces is rather sensitive to the values of the coupling constants and to the form of the ttNN 
vertex. Experimental uncertainties in the parameters of the TM vr — vr potential are assessed. They 
lead to uncertainties in theoretical predictions of the triton binding energy of about ±0.4 MeV. 
The low-energy limit of tt — o" and tt — uj potentials is performed. It defines the coupling constants 
of effective contact ttNNNN vertices, which are compared with the corresponding contact vertices 
of Chiral Perturbation Theory. 

PACS numbers: 21.30.-x, 21.30.Fe, 21.45. +v, 21.10.-k 
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I. INTRODUCTION 



Nucleon-nucleon (NN) and three-nucleon (3A^) potentials to be used in calculations of 
properties of few nucleon systems at low and intermediate energies are these days successfully 
modeled through meson exchanges. Modern A^A^ potentials based on this picture provide 
an impressive description (with x^/datum ~ 1) of the NN scattering data below the pion 
production threshold and of the deuteron properties . These potentials often employ 
the underlying meson-nucleon dynamics only in one-boson-exchange approximation. One 



simple formalism, but it was also sup- 
which showed numerous cancellations 



reason for this is clearly the desire for a relative^ 
ported by extensive studies of the Bonn group 
between classes of diagrams with two and more exchanged mesons. The meson-exchange 
approach with heavy meson (and nucleon) resonances has been supplemented recently by 
a model-independent construction from a low-energy realization of QCD: chiral perturba- 
tion theory (ChPT). While some theoretical uncertainties in the formulation of this theory 
persist jj], the construction of the A^A^ potential in this framework is already approaching 
maturity with the next-to-next-to leading order (NNLO) potential also describing 

the A^A^ data fairly well. Moreover, a recent paper ^ by the Jiilich-Bochum group provides 
a detailed comparison of ChPT with phenomenological potentials based on heavy meson 
exchange, and shows that the low-energy constants of ChPT appearing in the two-nucleon 
sector of the chiral Lagrangian agree reasonably well with those following from reducing the 
heavy meson operators to their contact form. Therefore a mutual positive influence between 
the two approaches should be expected for the future: the phenomenological potentials 
should provide a hint on the values of the chiral low-energy constants (which are not always 
easily obtained from the experiment) and the chiral constraints should be imposed on the 
phenomenological models at low energies. In this spirit some studies of ChPT 3A^ force 
were performed j^Q], although more extensive numerical calculations are still needed. To 
describe the data one might have to consider the chiral 3A^ force from the next order of the 

In the present paper we investigate the role of some short-range heavy meson exchanges 
in the 3A^ potential with respect to the binding energy of the 3A^ system. In addition to the 
dominant vr — tt component of the 3A^ potential Q,[l2|5 we include the force following from 
IT — p exchange and the terms involving n — u or ir — a exchanges from Ref. [lJ|. The 
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FIG. 1: Contributions to the 3A^ force. Diagrams for (a) nucleon Born contributions, (b) inter- 
mediate Roper resonance, and (c) vr re-scattering. Diagrams (a) and (b) involve one pion and 
one heavy- meson exchange. They are cahed "long-short range 3A^ forces". From diagrams (a), 
iterations of the A^A^ potential are removed, as described in Appendix Diagram (c) involves 
only pion exchanges. It is called a "long- long range 3A^ force". The meson momenta are defined 
as qi = p[- Pi, i = 1,2,3. 
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latter potentials were obtained both from diagrams with the Roper resonance A^*(1440) (Fig. 
^p) and from the nucleon Born terms (often imprecisely called pair- or Z-diagrams) where 
an off- mass-shell nucleon propagates in the intermediate state (Fig. ^). All contributions 
of the Born diagrams are of relativistic order. The parts with negative energy propagation, 
being genuine "Z-diagram" contributions, are included in the irreducible 3A^ potential. Also 
the part with off-shell nucleons propagating with positive energy cannot be omitted. These 
terms differ in general from the iteration of the A^A^ potential (which is not to be included 
in the 3A^ force and has to be subtracted from the total amplitude explicitly) and thus also 
contribute to the 3A^ potential. 

The purpose of this paper is to study the effects of the 3A^ forces based on pion-pion 
and pion-heavy-meson exchanges. It is discussed how strongly they depend on experimen- 
tal uncertainties in the determination of the phenomenological constants in terms of which 
these potentials are parameterized, and on the model dependence of the nucleon Born con- 
tributions. These inherent uncertainties should provide - for such semi-phenomenological 
meson exchange models - theoretical error bars of the calculated bound state energy of the 
trinucleon system. 

The latest version of the phenomenological Tucson-Melbourne (TM) potential 
TT — vr potential is parameterized in terms of three constants a, b, d. These constants are 
determined from vrA^ scattering data and we look at the variations of the triton binding en- 
ergy within the range of their experimental uncertainty, complementing thereby the usually 
studied sensitivity to the variation of the form factor cut-off parameter. 

The pion-heavy-meson exchanges contributions contain also some not very well deter- 
mined constants, in particular the coupling constants in vertices involving the Roper res- 
onance A^*. But besides that the Born terms also depend on the employed form of the 
ttA^A^ couplings. It is often assumed that the pseudovector (PV) coupling is preferable: 
this is based on the experience with modeling relativistic OPE A^A^ potentials where the 
pseudoscalar (PS) txNN coupling implies unrealistic enhancement of the intermediate neg- 
ative energy states. It is also claimed that a PV coupling is preferred because it does not 
require the non-minimal contact terms to conform with chiral symmetry in processes like 
pion electroproduction on the nucleon. 

However, in constructing the vrA^ scattering amplitude, or amplitudes of heavy meson 
production in pion absorption on the nucleon, one has to be more careful. In these cases one 
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cannot state that the use of PV couphng minimizes the contact Lagrangians: while it is true 
for the isospin-even ipil)^"^ vertex, the isospin-odd contact Weinberg- Tomozawa interaction 
is actually close to zero in the representation with PS coupling. If one tries to include the 
heavy mesons the situation becomes even less clear: to include them one has to extend the 
symmetry in a model dependent way, e.g. to require the local SU(2) x SU(2) symmetry 
dynamically realized (and broken) either in the Young-Mills fashion or using the concept of 
hidden symmetry realization. The construction of these (approximately) chiral symmetric 
Lagrangians including heavy mesons is thus model dependent. It does not allow to formulate 
the consistent chiral counting scheme and even at the tree level the dynamical content has 
not been sufficiently constrained by detailed analysis of wide range of hadronic processes, 
which they in principle should describe. 

In this paper we therefore consider only the Born amplitudes and the amplitudes with 
intermediate excitation of the Roper A^*(1440) resonance, which have been proposed in 
Ref. Q, and have never been included in realistic calculations of trinucleon bound states. 
Since we conclude from the reasoning above that there is no strong reason to prefer the 
Born contributions obtained with PV ttNN coupling, we calculate also their PS version. In 
contrast with Ref. [14] we include all these potentials in 3A^ Faddeev calculations exactly, 
i.e., without using perturbation theory. Their contributions to the binding energy are given 
individually. 

Finally, we relate some of these short-range 3A^ potentials to the corresponding coun- 
terterms from the chiral Lagrangian. We have deduced from the potentials with heavy 
meson exchange the effective contact low-energy four-nucleon-pion coupling constants and 
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15|. When 



attempted to relate them with the constants of NLO interactions of ChPT 
the heavy meson propagator is reduced to a point, it appears that the pion-heavy- meson 
3A^ potentials with PS ttNN coupling are closer to those obtained from the contact low- 
energy four-nucleon-pion NLO interactions of ChPT jlfll . Q| than their PV ttNN versions. 
Nevertheless, the comparison with chiral low energy constants is not straightforward, since 
some of the effective interactions we obtained by taking the point limit are not included 
in the NLO interactions of ChPT, and can be transformed to that form only after certain 
approximation. This approximation however does not seem to be numerically supported by 
results of our model calculations. 

This paper is organized into four sections. Section II contains an overview on 3A^ forces. 
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Section III presents the long-long and the long-short 3A^ forces and the numerical results. 
Section IV determines the low-energy constants and section V gives a summary and conclu- 
sions. 



II. BRIEF OVERVIEW OF 3N FORCES 

Chiral symmetry has been recognized as an important guideline for the construction 
of nuclear forces, into which the vrA^ amplitude enters as one of its building blocks, long 
before the advent of ChPT. The process of incorporating constraints from chiral symmetry 
breaking) into theoretical studies of the two-nucleon interaction was pioneered by G. Brown 



Iq in the early 1970's. Later, the same ideas were applied to the 3N force, for which two 
different approaches were developed in parallel: 

(i) One is based on the concept of partially conserved axial-vector current (PCAC) and 
the current algebra (CA) formalism, built from equal-time commutation relations for vector 
and axial-vector currents. It underlies the well-known Tucson-Melbourne (TM) two-pion- 



exchange 3A^ force ll|, represented in Fig. It goes beyond the static P-wave A 
contribution to the vrA^ amplitude, the only building block considered in its predecessor, 
the Fujita-Myazawa 3A^ force. Current algebra and PCAC provide an elastic pion-nucleon 
scattering amplitude which includes the pion-nucleon cr-term, a direct measure of chiral 
symmetry breaking, which can be extracted from experiment. 

(ii) The other approach stems from the so-called "effective chiral Lagrangians" for the ttN 
system. The first example of this approach is the Gell-Mann and Levy linear sigma model. 
Another example is the chiral Lagrangian with pseudo-vector vrA^ coupling, which underlies 



the construction of the so-called Brazilian 3A^ force in Ref. 117 1 . More recently, ChPT, based 
on the Weinberg Lagrangian |18l| supplemented bv multinucleon contact terms, became a 
systematic way to approach two and 3A^ forces and was applied to the description of 
low-energy hadronic physics. 

In particular, it was realized in Ref. that although the construction of the TM two- 
pion exchange 3A^ force employs chiral constraints for ttN scattering through CA, it missed 
further constraints, which arise from embedding that amplitude in the 3A^ system. The 
conclusion was that the TM form of the two-pion-exchange force contained a spurious term 
corresponding to contact terms between two nucleons and pions {NNtttc). This spurious 
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term can be generated also in ChPT with the help of a pion field redefinition. But at the 
same time an additional two-nucleon contact term arises {NNNNtc), which has not been 
taken into account by the TM group and which exactly cancels the first spurious contact 
{NNtttt) contribution. As a net result, the so-called c term of the TM force should be 
dropped and the so-called a term is modified. We present here the effects of the chirally 
imposed changes to the TM force on the triton binding energy. A previous calculation of 
these effects j2oj used a less realistic A^A^ interaction and a variational numerical method. 

The CA program was extended to describe also nN-pN transitions, from which a tt — p 
exchange 3A^ force can be constructed. It seems natural to include such a mechanism, given 
the important interplay between vr and p exchange in two-nucleon potentials. To model the 
TT — p 3N force, one can use vector-meson dominance to access the p-analogue of the off- 
mass-shell pion electroproduction. Chiral symmetry is in this case supplemented by gauge 
invariance to constrain the Ward identity amplitude. 

The interest in the short-range 3A^ forces increased in mid 1990's, when it was found 
that the lack of binding energy of the triton is not the only experimental signature of a 3A^ 
force. According to 21j a spin-orbit structure of the 3A^ force, not present in the standard 
two-pion exchange based forces, may solve systematic deviations between the theoretical 
and the experimental maximum of the nucleon vector analyzing power Ay (6) in elastic nd 
scattering below 30 MeV. The TM or the Brazilian 3A^ force either produces no sizable effect 
or even worsens the disagreement with the data. But the study of Hiiber et al. j^l indicated 
that the short-long range forces with a structure similar to that given by vr — a and -k — uj 
exchanges with an intermediate Roper resonance A^* and suggested by the NLO ChPT 
jisl ] could potentially improve the description of Ay. Similarly, a purely phenomenological 
spin-orbit 3A^ force yields a noticeable improvement of the description of this observable 

More recently, Meissner et al. j3| actually claimed that the Ay puzzle can be resolved 
by a chiral A^A^ force alone. However, the A^A^ interaction employed does not match the 
high quality of contemporary descriptions of A^A^ scattering data, in particular, for the 



crucial triplet P-waves 
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21[. The NLO chiral potential gives a qualitative description of 



the phase shifts only for very small energies (< lOMeV), but for a truly quantitative fit one 
has to include NNLO corrections . Entem and Machleidt showed that no low-x^ A^A^ 
potential, neither one based on heavy-meson exchanges nor one constructed from ChPT, 
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can solve the Ay puzzle, since a good fit to the A^A^ phase shifts (in particular in the ^Pj 
channel) does not allow the variation needed to fix this spin observable. In a recent paper 
of Epelbaum et al. j3| it was found that NNLO A^A^ and leading order 3A^ chiral potentials 
still do not solve the Ay problem completely. Also the Urbana 3A^ force (a Fujita-Myazawa 
force with a short-range regularization determined phenomenologically through a fit to the 
triton binding en ergy ) does not provide a solution. An interesting alternative was suggested 
by Canton et al. j2J], who constructed an additional 3A^ force by reducing tt — 3A^ equations 
below the pion threshold. Their force is very similar in spin-isopsin structure to the vr — a 
and vr — force with intermediate Roper resonance discussed above, but it is multiplied by 
a part of an energy- dependent A^A^ T-matrix. 

In addition to the mentioned nd scattering problem, a good description of neutron-rich 
light nuclei is only possible 2^ when a phenomenological 3A^ interaction with three pion 
exchanges is added to the original Urbana force 2^. Finally, let us mention that while 
realistic A^A^ potentials alone predict an equilibrium density for symmetric nuclear matter 
that is too high, phenomenological 3A^ forces can change it in the right direction j^]. 

This overall picture indicates that the long-range tt-tt 3A^ force needs to be complemented 
by other physical mechanisms. 



III. LONG-LONG AND LONG-SHORT RANGE 'iN FORCES 

The phenomenology of the A^A^ interaction shows that the most important meson ex- 
changes in every realistic OBEP are exchanges of the pion, of the vector p and uj mesons, 
and of the (fictitious) scalar a meson. Therefore, one should also investigate the role of such 
exchanges in the modeling of a 3A^ force. 



A. vr — TT exchange force 

Given the importance of the one-pion-exchange potential (OPEP) in determining the 
long-range part of the A^A^ interaction, it was natural that the first 3A^ force considered 
was the two-pion exchange. The corresponding diagram is shown in Fig. where the 
blob represents all possible pion re-scattering processes. The iteration of the OPEP has to 
be subtracted from this Feynman diagram, since it is generated automatically by iterating 
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the NN potential in the Faddeev equation. The ttN amphtude tT^N entering the 3A^ force 
involves virtual pions. This off-mass-shell continuation is constrained by the soft pion low- 
energy theorems. In this paper, we adopt the (modified) TM version of vr — tt force, which 
has been constructed to explicitly conform with these requirements. 

The TM vr — vr 3A^ force was generated by applying the Ward identities of CA to 

the amplitude of axial-vector current scattering on the nucleon. Using PCAC, the resulting 
TfN amplitude, in the even and odd isospin channels, was written as an expansion in powers 
of the momentum transfer t and the crossing- variable u = {s — u)/4M. Then the "reality 
test" was applied with the successful result that CA predicts the first four coefficients of 
the expansion in good agreement with the empirical subthreshold expansion coefficients 
obtained from dispersion relations. Therefore, it was possible to use accurate empirical vrA^ 
data to construct the nearly model-independent TM force. The dispersion analysis has been 
updated over time with the inclusion of new data, therefore the coefficients of the TM force 
have also evolved. 

In terms of physical mechanisms, the non-Born part of the TM vrA^ amplitude contains 
terms that can be interpreted as t-channel a and p exchanges, as well as contributions from 
the A resonance in the intermediate state. The contributions of the subtracted nucleon Born 



diagrams were a naly zed in detail in 
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29| and it was later shown numerically that they are 
relatively small 3Q], at least as far as their contribution to the triton binding is concerned. 

As argued in ;15], the TM ttN amplitude, although constructed from CA and obeying the 
chiral constraints in the nN sector, should not be simply attached to two additional nucleon 
lines, since that would be inconsistent with the chiral counting for the 3A^ potential. More 
precisely, the TM vrA^ amplitude can be derived from the usual chiral Lagrangian after the 
pion field re-definition 

7z' = 7r{l- c)N^N , (1) 

where c ~ cr/^rn^f^) is the constant of the so-called TM c-term. The original chiral La- 
grangian in the "natural" representation acquires after such redefinition two additional terms 
(relevant up to the order considered) 

A£(^) = -cN^N{{7z'd^d''7r' + ml7r"') + ^[{VN^)-(TT-n'N + NW-7z'(T-iVN)]}. (2) 

The TM amplitude includes the term corresponding to the first part of AC^^^ (the c-term 
proportional to + q^), but in the construction of the TM 3A^ potential the second term 



(which involves two nucleons) was not considered. Therefore, according to Ref. |15|, the TM 
c-term should be dropped and its a-term should be replaced by a' = a — 2m^ c. The TM 
71 — 7T force modified in this way was recently adopted in [2^ and labeled TM'(99). However, 
we would like to point out that it is not correct to transform away the whole c term, because 
it includes the contribution from the Born diagrams- that contribution arises naturally also 
in ChPT in its usual representation. Nevertheless, since this residual c-term is rather small 
and the tt — vr force is not the main subject of interest of this paper, we avoid introducing 
further redefinitions of the TM force and we adopt TM'(99) for our numerical calculations. 
Thus, the vr — vr potential used in this paper reads 

-.2 



W^l) = -(27r) 
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4^3 



(0-2 ■ qa) (0-3 ■ qg) A^(q2)A^(q3) 



X {(t2 ■ T3) [a' - 6q2 ■ qs] - (m ■ T2 X T3)d {i(Ti ■ qa X q3)} + 2^3, (3) 



where we use for the meson momenta qj = — Pi, so that they are always pointing 
away from the "active" nucleon as shown in Fig. ^ (to avoid confusion, we note that the 
different convention with the first pion incoming and the second outgoing from the active line 
rescattering "blob" is also often used in discussion of this potential). The overall momentum 
conservation is then given by qi + qa + qs = 0. The functions AB(q) are defined to contain 
the propagator function of the meson B and the square of the strong BNN form factor 
■^-B(q) = ^BNNi^^)/^^^ + ^^)- The full potential contains additional four terms following 
from (jS)) by cyclic permutations. 

The parameters of the TM'(99) force are given by 

,2 



cr 



- 2m! 
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b 
d 



2 



m 



F+(0,m^H-)-- 



Arn^ 
a 

J TT 



+ Knn{0) 



a 

f TT > 



B-{0,0) ^ g 



2m 



Am? 



(4) 
(5) 
(6) 



Here, a is the pion- nucleon sigma term, g is the ttNN coupling constant, the pion 
decay constant, F'^{i^, t) and B^{u, t) are isospin-even non-spin flip and isospin-odd spin-flip 
t-channel vrA^ amplitudes, respectively, with the nucleon pole term subtracted Note 
that we distinguish between the mass of the charged pions, mT^+ = 139.6 MeV, and the 



isospin averaged pion mass = 138.0 MeV. Although their difference is small, it affects 



10 



^iTNN/'mn 



4.1 



-1.203 



5.0 



-1.154 



5.8 



-1.127 



6.5 



-1.112 



7.1 



-1.101 



TABLE I: Expansion coefficient a' of the vrN amplitude used in the TM vr — vr force as a function 
of the ttNN vertex cut-off parameter Ayr at at, in units of the isospin averaged pion mass. The other 
coefficients b and d do not depend on AttAtat and have the values m^b = —2.801 and m^d = —0.754. 

the 3N force parameters in a noticeable way. The vertex form factor 



depends on the cut-off parameter A^^nn, which in the original TM force was taken to be 5.8 
in order to be consistent with a Goldberger-Treiman discrepancy of 3% determined at 
the time. Based on more recent data, the Goldberger-Treiman discrepancy shrank to about 
2%, which corresponds to A^A^Af close to 7.1 m^. On the other hand, A^A^Af is frequently 
tuned to reproduce the triton binding energy in calculations that include only the vr — tt 
exchange part of the 3A^ force. Such calculations require much lower values, close to 4 m^r. 
Clearly, the functional form of the form factor is not much constrained by the knowledge 
of the coupling constant g FT^^N^'f) at the two points = and = and thus does 
not allow a strict determination of a cut-off mass. We keep therefore the by now traditional 
reference value A^A^Af = 5.8 m.^ as our standard one, but also vary it then between the limits 
indicated above to study the sensitivity of the results. In particular, we investigate the 
question whether the short-range forces are able to reduce the strong cut-off dependence, as 
one may expect from arguments of ChPT. 



Adopting the values ^ g"^ = 172.1, F+{0,ml+) = -0.05 m;^, 5^(0, 0) = 8.6 m;^, 
-j2 = 1-40 and ATrAfA^ = 5.8 m.,^, we obtain the coefficients given in Table HI Note that 

J TV 

a' depends weakly on Attatat, a dependence often ignored in practical calculations that vary 
the cut-off mass. 




(7) 
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While the rather strong dependence of the 3A^ binding energy on A^^nn has been recog- 
nized as a source of significant uncertainty, httle attention has been paid to uncertainties 
originating from the experimental errors in the other parameters of the TM vr — vr exchange 
force. Since the first publication of the TM force, the experimentally determined values of 
g"^, F+(0,m^+), 5^(0,0), and have changed several times, leading each time to updated 
force parameters a', b, and d, according to Eqs. (H)) to ©• We estimate here how these 
experimental uncertainties propagate into uncertainties of the triton binding energy Ef. 
For simplicity, we introduce the dimensionless variables 



y = g\ z = F+(0,m2 + )m,+ , n = S-(0,0)m2 + 
A 



A 



as well as the dimensionless force parameters 



, 9 = — , 

m 



(8) 
(9) 



a = a! fJL , h = hjjL^ , d = djj? 



(10) 



Using this notation, Eqs. (jlj) through (jH)) become 

A^-r 



a = xr \ 1 — 2r — 2- 



A4 



1 3 



b = -2r'^{x - z) , 

J 1 3 1 2 
d = -^q y--r qu. 



(11) 
(12) 
(13) 



Adding (independent) errors in quadrature, and taking into account the relations ((TT} 
through ()13p . we arrive at the following expression for the square of the uncertainty in 
due to the uncertainties Ax, Ay, Az, and Am in x, y, z, and u: 



(AEt) 
"1 



r 1 - 2r^ - 2 



A2 - 1\ dE, 



A4 



'dE, IdEt 



da 2 dd 



da 



2r 



,dE, 



1 2 



db 



iAxf+ 



>dE^ 
db 



1 2 



{Azf + 



1 2 dEt 



1 2 



(An) 



(14) 



Note that it is simpler to vary the three constants a, b, and J instead of the four experimental 
values X, y, z, and u, but one has to keep in mind that the former are not independent from 
each other. 
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Reid Paris Nijmegen 93 Bonn B 

dEt/da 0.080 0.090 0.095 0.095 

dEt/dh 0.725 0.745 0.845 0.750 

dEt/dd 0.545 0.470 0.540 0.505 

/^Et 0.377 0.389 0.440 0.393 

A£;fPP 0.362 0.372 0.423 0.375 



TABLE II: Numerical results for the partial derivatives of the triton binding energy with respect 
to the TM TT — vr force parameters and resulting uncertainties according to Eqs. (fTl)) and ((T3|l . 

The partial derivatives dEt/da, dEt/dh, and dEt/dd, are calculated numerically for the 
standard set of parameters of Table U at A = 5.8. As estimates of the experimental uncer- 
tainties in X and z we use the values given in Ref. [2^, Ax = 0.25 and /\z = 0.05. As a 
reasonable estimate for Ay we choose the difference between the current value y = 172.1 
and the one used in the original TM force, y = 179.7, yielding Ay = 7.6. The value u = 8.1 
is given in [2^ without indication of the error. Assuming that the specified digits are indeed 
significant, we set Au = 0.1. 

The numerical results for the partial derivatives of Et with respect to a, b, and d, and the 
corresponding values of AEt are shown in Table [H] for the TM vr — vr force in combination 
with various A^A^ potentials. In each case, AEt is roughly 0.4 MeV. This is a significant 
value which clearly shows that the dependence on Attatat is not the only source of uncertainty 
in predictions of Et. 

A closer inspection of ()14|) reveals that AEt is almost completely dominated by Ax, the 
uncertainty in the pion-nucleon sigma term. One can easily derive the approximate form 

d F 

AEr^2^Ax, (15) 

the results of which are also shown in Table [H] and come very close to the full result, 
confirming that the uncertainties in the other parameters are secondary as long as the pion- 
nucleon sigma term is not determined with significantly better accuracy. 

Since fjl5p is a lower limit of (fT^ one also has to conclude that currently predictions of 
the triton binding energy employing the TM force cannot be made with a better accuracy 
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than plus or minus 0.4 MeV. Clearly, adding other contributions to the 3A^ force will further 
increase this uncertainty. 



B. TT — p exchange force 

In the A^A^ potential, p-meson exchange provides important contributions to the tensor 
and spin-orbit components. Its role is enhanced by the large anomalous pNN coupling Kp, 
connected (via VMD) to the anomalous isovector magnetic moment of the nucleon. The 
71 — p exchan ge w as therefore always considered the next most important part of the 3A^ 



force 
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29l |30| after vr — tt exchange. In this paper we will take another look at its 



contributions stemming from the nucleon Born diagrams. The TM vr — p exchange force 
also includes processes with intermediate A resonance excitations, which numerically turn 
out to be equally important. In our calculations the tt — p potentials with intermediate A 
resonance are included in the form (and with the corresponding parameters) specified in 
Ref. 0. 

The nucleon Born contributions to the vr — p 3 force are derived in Appendix El They 
depend significantly on the form of the nNN coupling. For the PV nNN coupling one gets 



W. 



PV+I 
■KpT ' 



1^ 



■{271 



Qi] 



X 



0-2 • q2 A^(q2)Ap(q3) +2^3 



W. 



PV+/ 
npS ^ 



0, 



(16) 
(17) 



whereas for the PS nNN coupling 

4 

2^2 



-(2^)^ 



(r2 ■ T3)cri ■ qs 0-2 ■ q2A^(q2)Ap(q3) +2^3, (18) 



PS-/ 



.(2^.3^ 9p 



Arrfi 

ITi ■ T2 X Tg) [cTi ■ Q3 + (1 + Kp)i(Ti X CTg ■ qg 



4m^ 

(72 ■ q2 A^(q2)Ap(q3) + 2^3, 



(19) 



where again qj = — Pi, Qi = Pi + Pi- The subscripts S and T refer to the exchange 
of the space and time components of the p field, respectively. For brevity, we will call the 
corresponding terms "spacelike" and "timelike". The superscripts PS and PV correspond to 
the type of 71NN coupling, the superscripts "+" stand for the contribution of the positive 
energy nucleon Born diagrams as derived in Appendix \K\ superscripts "-" denote the "true" 
pair (or Z-diagram) terms. 
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Meson mass (MeV) ^ ^ cut-off 

vr 138 13.6953 2.1664 A^a^tv = A^Ar7V.=5.8m^ 

p 768.3 0.81 ApNN (Dirac)=12m^ 

ApjvAr(Pauli)=7.4m^ 

a 584 10.3251 0.6 A^^tvat = A^attv. =1995 MeV 

LO 782.6 24.5 1.4237 A^atat = A<^ArAr.=1850 MeV 



TABLE III: Masses, couplings, and cut-ofF parameters that appear in the 3A'^ forces. For the p and 
the oj mesons there are also the tensor /vector coupling ratios Kp = 6.6 and = 0.0, respectively. 
The parameters of the vr and p mesons are taken from the TM'(99) force and from Ref. 3^, the 
parameters of the uj meson from the Bonn B potential. The a meson parameters are determined 
in Appendix IbI 

In the numerical calculations we keep the tc — p force parameters of with the 

exception of which is updated to the value of TM'(99) (see Table ITTT]) . 

The terms proportional to Qj are non-local and have been omitted so far from 3A^ poten- 
tials, mainly because of the difficulties associated with performing calculations with non-local 
interactions in coordinate space. Working in momentum space, we are in principle not ham- 
pered by non-localities. However, in this paper we focus on the relation between traditional 
meson-theoretical 3A^ forces and the 3A^ forces derived from ChPT where non-local terms 
are discarded. Hence it would not be useful to keep those terms in our calculations, and we 
neglect them as well. 

The above results illustrate one important point: One cannot rely on the nucleon Born 
terms alone to yield the most important contribution to the 3A^ force. In particular, it 
would be incorrect to claim that since the PV ttNN coupling is "more consistent" with 
the important chiral symmetry requirement, one should use the above potentials VTjpj^ 
and W^^g^. Indeed, the construction of the it — p 3N force by the TM group fl3^ via the 
extension of CA to include vector meson dominance, has shown that its most important 
part - the so called Kroll-Ruderman (KR) term - follows (in the theory with PV coupling) 
not from the nucleon Born term, but from an additional contact term required by gauge 
invariance of the pion photoproduction amplitude. This is, of course, very similar to the 
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situation for the leading order isovector meson exchange currents: there the nucleon Born 
term with PS couphng gives the most important contribution (which is supplemented by 
the pion-in-fiight diagram) , whereas the Born term for P V coupling does not contribute and 
gauge invariance requires the presence of the contact interaction. The KR contact term 
happens to be opposite in sign to the chiral contact term, arising from the chiral rotation 
from the PV to PS coupling and introduced in Appendix |X1 Thus, it appears reasonable to 
adopt for the tt — pS exchange potential just that KR term, i.e., to take 

w^fpf = -w:;f = w:^s (20) 

But then it would seem natural to take for the timelike exchange the corresponding W-j^pT 



W^pj^, which is enhanced by the large factor 1 + Kp compared to W^Jj^ . Note, that in 14 1 
(in which the Born contributions to the tt — pS where first derived) an expression identical 
to W^Jpj^ is listed for the timelike part, while W^p^ (the KR term) is given for the spacelike 
part ^. 

We do not intend to imply by the discussion above that the PS tcNN coupling is to be 
preferred in this context. In fact, trying to make any preference for the Lagrangian involving 
heavy mesons is meaningless (unless one attempts to extend the global chiral symmetry to a 



local one and builds a model which provides all relevant contact terms 



3l|). We only want 



;o point out that more simplistic approaches (like the one of this paper and that of Ref. 



14| : taking only the Born terms with a hope that they contain the most important effects) 
do not give an unambiguous prediction. Besides, the p — T contribution can be comparable 
to others when it is enhanced by the large factor 1 + /tp, as it happens for PS coupling. 



C. TT — a and tt — uj exchange force 

Since a and u meson exchanges play an important role in OBE models of the NN 
interaction, they contribute naturally also to the 3A^ force. Short-range tt — a and tt — u 3N 
potentials were introduced in Ref. jl^- They are derived from diagrams with intermediate 
positive-energy nucleons with PV ttNN coupling. 

^ There are, unfortunately, two misprints in Eq. (2.8) of Ref. a factor of (1 + Hp) is missing in the 

second term and the first term has an extra factor 2. 
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Unlike in the vr — vr force, the Born terms in vr — cr and vr — potentials (Fig. are 
rather large. As the tt — p potentials from the previous section, they depend on the type 
of ttNN coupling. We also include additional a and u exchange contributions generated by 



excitations of the intermediate nucleon to the Roper resonance (Fig. H^)- 

The TT — a and n — u potentials corresponding to the nucleon Born diagrams follow 
immediately from the expressions developed in Appendix ^ For vr — a exchange we get 

</+(!) = +(27r)30(Ti ■ T2) a, ■ qs ^2 ■ q2 A.(q2)A.(q3) +2^3, (21) 

2 2 

W^^-il) = -(27r)3|||(ri • T2) a, q.cr,- q^ A.(q2)A.(q3) +2^3, (22) 

and for vr — exchange 

iyPJ+(l) = -(27r)30(ri-T2)^i-q3a2-q2A.(q2)A^(q3) + 2^3, (23) 
<J-(1) = -(27r)^ ^'^'^^+ "^"^ (n • T2) a, ■ qs ■ q^ A.(q2)AUq3) + 2^3. (24) 

In Ref. 



14| . only the potentials for PV coupling are given. Since k,^ is very small, the 
results for it — u do not depend much on the type of the nNN coupling. However, the 
71 — a potentials do differ. For PV coupling, the vr — cr and n — lj potentials W^^~^{1) and 
W^^~^{1) have identical structure but opposite sign. Therefore, as in the case of the NN 
interaction, a strong cancellation between the vr — a and n — u 3N interactions occurs. For 
small momenta, the sum of these potentials is proportional to 

4-4' (25) 

just as for the corresponding A^A^ potentials. For the potentials with PS coupling, Eqs. 

such a cancellation does not take place. However, they would cancel if the potential 



did not change much when the momentum is replaced by the — q2- As discussed 
later, this is exactly the momentum replacement needed to extract the low-energy constants 
(LECs) of ChPT. Therefore we calculate the tt — a potential in both forms. The extent of 
the cancellations between a and uj exchanges depends also on the numerical values of the 
masses and coupling constants of the a and uo mesons. 

While we are trying to parameterize the 3A^ forces as consistent as possible with the A^A^ 
potentials they are combined with, we face a problem with the a meson. The only A^A^ 
potential with explicit a exchange we are using is Bonn B. Its cr meson, however, is not a 
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pure scalar-isoscalar particle, but has also a scalar-isovector component. This is reflected in 
the fact that the a mass, the aNN coupling constant and cut-off mass is different in A^A^ 
isospin and 1 channels. 

In order to get a true OBE representation of a-exchange, we keep for the 3A^ potential 
only the part of the Bonn a meson which corresponds to pure isoscalar exchange. The 
parameters of our a are then related to the Bonn a such that a simple OBE form, consistent 
with the way it is implemented in the 3A^ forces, is obtained. Details of this procedure can 
be found in Appendix iBl 

The forces originated from excitations of the Roper resonance read 

W:AI) = -(2^)' 2(J-'^^]».2 (^1 • ^1 • q2 ^2 • q2 A.(q2)A.(q3) +2^3, (26) 

W:M) = +(27r) % ^f^'^^*- ^ in ■ T,) a, ■ q2 0-2 ■ A.(q2)A^(q3) + 2^3. (27) 
2(m — m)m'^ 

Here m* is the mass of the A^*(1440) resonance. We would like to stress that the form of 
these potentials depends on the structure of the NN*a and NN*uj vertices. We adopted 
the very simple choice used in Ref. Q]. The form of these effective vertices follows from 
the quantum numbers of the baryon state A^* and of the a and uj mesons, neglecting any 
possible dependence on the substructure of these particles. Even for the simple choice of the 
A^* vertices of , theoretical predictions of the corresponding coupling constants (and their 
ratios) from quark models are hampered by the uncertain quark content of A^*, whereas their 
extraction from experimental data can be done only in an indirect and model-dependent way. 

The strength of the uNN* cou plin g is calculated from the partial decay width of A^* 
N + -K. The recent determination based on the value r(A^* Nix) = 228 MeV gives^ 

ZzEpi = 0.0117, / = /^^^. = 5.22 . (28) 

The couplings for the aNN* and ujNN* vertices are more difficult to pin down. The 



coupling constant of aNN*, extracted from the partial decay width of A^* — >■ A^+ {'n"i^)s-wavej 



1=0 
S—w 

The "(T meson" 



depends critically on the assumed mass and width of the a resonance 
(with zero width) employed in the parameterization of the OBE potentials simulates not only 
{'^''^Ys^wavey but also Other scalar-isoscalar exchanges. Such a phenomenological a exchange 



^ Equation (3.3.) of Ref. Q| is incorrect and the corresponding value of is overestimated by a factor 

of 3. 
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was used in a recent analysis 



33| of the inelastic scattering a + p — > a + X, from where the 



effective coupling constant of aNN* was extracted as 

1.33, 



9^ 



(29) 



with m„ = 550. This is much larger than typical values obtained from the Roper resonance 
decay |32l |. 

On the other hand, Ref. 



of Ref. 



3- 



extracts 



47r 



0.1, which is even smaller than all values 



Given this wide spread of coupling constants in the literature, we adopt the 



intermediate = 0.6 as our standard value, but calculate also the two extreme cases in 
order to observe the sensitivity of the results with respect to the choice of this parameter. 
For uNN* we follow Ref. 



14| and determine from the ratio 



9u> = 9u 



9a 



(30) 



which follows from a naive constituent quark model. The sum of the 3N potential with 
intermediate Roper resonance and it — a and ir — u exchanges is then roughly proportional 
to 



9a 



9l 



9: 



. . , • (31) 

Therefore we can expect the same amount of cancellation between a and uo terms in 3A^ 
forces with Roper excitations as in NN potential and 3A^ Born contributions. However, 
since the simple scaling rule of coupling constants is theoretically not very well founded, we 
also use sets of parameters that do not satisfy (jHUjl . 



D. Numerical results 



We calculated the triton binding energies and 3A^ wave functions for Hamiltonians con- 
taining both NN and 3A^ potentials by solving the non-relativistic 3A^ Faddeev equations 
in momentum space exactly. 

The 3A^ potentials are expressed in terms of Jacobi momenta in the cm frame of the 3iV 
system, and decomposed into partial waves in a basis of jj coupling states. Because of the 
complexity of the resulting numerical problem, we restrict ourselves to those partial waves 
in which the NN pair total angular momentum does not exceed 2, which corresponds to 
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Reid Paris Nijmegen 93 Bonn B 

3NF Et AEt Et AEt Et AEt Et AEt 



no 3NF 


-7.230 






+7r7r(a') 


-7.279 


-0 


049 


+7r7r(6) 


-8.739 


-1 


460 


+7r7r((i) 


-9.100 


-0 


361 


+TTp{KR) 


-9.017 





083 


+vr/3(A+) 


-8.849 





168 


+vr/)(A-) 


-8.747 





102 


+7Tp{T) 


-8.772 


-0 


025 




-8.273 





499 


+TTa{N*) 


-8.711 


-0 


438 


+-K0J{Z) 


-9.213 


-0 


502 


+^uj{N*) 


-8.735 





478 



-7.383 




-7.756 


-7.439 


-0.056 


-7.811 


-8.939 


-1.500 


-9.471 


-9.220 


-0.281 


-9.782 


-9.118 


0.102 


-9.635 


-8.961 


0.157 


-9.464 


-8.821 


0.140 


-9.285 


-8.850 


-0.029 


-9.316 


-8.213 


0.637 


-8.663 


-8.610 


-0.397 


-9.145 


-9.380 


-0.770 


-9.977 


-8.898 


0.482 


-9.370 





-8.100 




-0.055 


-8.159 


-0.059 


-1.660 


-9.624 


-1.465 


-0.311 


-9.847 


-0.223 


0.147 


-9.672 


0.175 


0.171 


-9.506 


0.166 


0.179 


-9.325 


0.181 


-0.031 


-9.352 


-0.027 


0.653 


-8.658 


0.694 


-0.482 


-9.055 


-0.397 


-0.832 


-9.956 


-0.901 


0.607 


-9.524 


0.432 



TABLE IV: Triton binding energies and their differences (in MeV) calculated for various model 
Hamiltonians with different A^A^ potentials and contributions to the 3A^ force added consecutively. 
All ttNN vertices in the 3A^ forces of this table are calculated in PV coupling. The columns 
labeled Et show the triton binding energies, while the ones labeled AEt indicate the differences 
between the binding energies of consecutive rows, indicating the effect of the corresponding 3A^ 
force component. 



18 different three-body channels. Details of the formalism and the numerical methods are 
described in Ref. 3^ . 

The main purpose of these calculations is to compare the effects of the various contribu- 
tions to the 3A^ force which were described in the previous sections on the triton binding 
energy. In order to study also the dependence of our results on the A^A^ interaction included 



in the Hamiltonian, we used a number of different A^A^ potentials, namely Reid j38|], Paris 
0, Nijmegen 93 Q, and Bonn B 40|. 



In Tables HVl and fVl we show the triton binding energies for various Hamiltonians where 
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the 3A^ forces of the previous sections are added successively. The binding energy differences 
depend somewhat on the order in which the potentials are added, therefore we later also 
show their expectation values. 





Reid 




Paris 




Nijme^ 


fen 93 


Bonn B 




3NF 


Et 


AEt 


Et 


AEt 


Et 


AEt 


Et 


AEt 




-8.859 


-0.112 


-8.953 - 


-0.132 


-9.425 


-0.140 


-9.452 


-0.127 


+vra(Z) 


-10.492 


-1.633 


-10.544 - 


-1.591 


-11.431 


-2.006 


-10.879" - 


-1.427 


+tt(t[N*) 


-11.219 


-0.727 


-11.264 - 


-0.720 


-12.215°^ 


-0.784 


-11.600" - 


-0.721 


+ TTUJ{Z) 


-13.680 


-2.461 


-15.367 - 


-4.103 


-20.952 


-8.737 


-15.811" - 


-4.211 


+7rw(iV*) 


-12.674 


1.006 


-14.134 


1.233 


-16.869 


4.083 


-18.345 


-2.534 



°In addition to the specified result, another unphysical, deeply bound solution was obtained in this case. 

TABLE V: Triton binding energies and their differences (in MeV), as in Table ITVl but with ttNN 
PS coupling in the 3A^ forces. The binding energy differences in the first row are calculated with 
respect to the corresponding entries in the row labeled +TTp{A-) of Table HVl Note that only 
for irpiT) and Tra{Z) the 3A^ potentials using PS and PV coupling actually differ. The effects of 
the other, unchanged 3A^ potentials on the binding energy are amplified compared to the case of 
PV coupling, since the 3A^ wave functions are significantly altered by the 7r/>(T) and TTa{Z) PS 
potentials. 



First we make a few comments to the results with the n — n 3N potential only. In our 
calculations ^30\ with the old version of the TM force TM(93), which contains the c-term and 
uses somewhat different values of the constants a, b and d, we obtained ii^t(Reid) = —8.904 
MeV and i?t(Paris) = —9.060 MeV. This means that transforming the c-term away and 
using the new version TM'(99), as in Table HVi brings an extra binding of about 200 keV. 
The first line of Table IVII shows that, even after removing the singular c-term, the TM 
TT — TT force remains highly non-perturbative. Also the cut-off dependence of Et (see Fig. |2)) 



remains about the same (see Refs. |20|, |30|) as with the old version TM(93 i. 

Most of the new short-range potentials considered in this paper contribute individually 
as much as (or more than) the TM it — p short-range forces considered before. The only 
exception is the p — T term, also not considered before. It is the only attractive part of the 
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3A^ potential 



no 3A^ force with vrvr with vrvr+vrp with TnT+irp+ira+TTUJ 



vrvr 


-1.837 


-1.281 


-2.572 


-2.190 


-2.274 


Tip 


0.370 


0.139 


0.450 


0.293 


0.317 


ira 


0.240 


0.275 


0.748 


0.605 


0.881 


TTUJ 


-0.288 


-0.228 


-0.758 


-0.594 


-0.968 



TABLE VI: Matrix elements of the 3A^ potentials of the first column calculated with different wave 
functions. For comparison, the column labeled AEt shows the binding energy differences according 
to Table ITVl All wave functions are calculated with the Paris potential and various 3A^ forces in 
the Hamiltonian. 

71 — p 3N force, but it is about 5 times smaller than the KR term. Only in PS coupling the 
p — T term is enhanced by the factor 1 + Hp and becomes comparable to the KR. 

The a and u exchange contributions are rather large. The Born terms are repulsive in 
case of the a and attractive for the uj, but when PV ttNN coupling is used they cancel in 
a similar way as the corresponding A^A^ potentials. For the PS coupling, the it — a Born 
contribution turns attractive and together the ir — a and n — uj exchanges lead to strong 
overbinding, as Table |V| demonstrates. Nevertheless, it would be premature to rule out the 
PS coupling solely based on these results: we would like to remind that other non-Born 
contact terms have not been investigated so far and are likely to be as important as the 
Born term considered here. 

One example of such contributions are forces with an intermediate Roper resonance. We 
observe that individually they are also rather large. However, as we have already pointed 
out above, there are significant uncertainties in the coupling constants involving A^*, and 
the effect of those potentials can vary accordingly. Nevertheless, it is remarkable that even 
for our smallest value of g* the tt — aN* term is comparable to the n — p A forces, as can be 
seen by comparing Tables HVl and Ivnl Although the effect of the tt — aN* force scales with 
g*, the overall contribution of vr — aN* and vr — uN* forces does not change much, since 
we assume the constituent quark model ratio g*/g^, which leads to a strong cancellation 
between them. Removing this constraint, as for instance in the last line of Table IVTTl can 
have an effect on the overall contribution of vr — cr and n — oo forces comparable to the 
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11 = 0.1 fi = 0.6 i = l-33 



3NF 




AEt 


Et 


AEt 


Et 


AEt 


vrvr + VT/O 


-8.850 




-8.850 




-8.850 




+7ra(Z) 


-8.213 


0.637 


-8.213 


0.637 


-8.213 


0.637 


+7ro-(iV*) 


-8.369 


-0.156 


-8.610 


-0.397 


-8.822 


-0.609 


+7rw(Z) 


-9.073 


-0.704 


-9.380 


-0.770 


-9.652 


-0.830 


with 5* = 5^ g 














+7ru;(A^*) 


-8.883 


0.190 


-8.898 


0.482 


-8.910 


0.742 


with |i = 0.24 














+7rw(7V*) 


-8.882 


0.191 


-9.174 


0.206 


-9.431 


0.221 



TABLE VII: Triton binding energies and their differences (in MeV), as in Table HVl for the Paris 
NN potential and 3A^ potentials with different coupling constants (7* and 5*. The first two lines 
are taken from Table Hvl to define reference values. The contribution of TTa{N*) increases with 
increasing g*, and so does the contribution of ttuj{N*) in the fifth line for a fixed ratio g^/g%, 
keeping the binding energy almost constant. When ^r* is kept fixed (last line), these cancellations 
no longer take place and the binding energy increases. 

transition from PV to PS coupling in Z-digrams: the cancellation between vr — aN* and 
TT — ujN* forces could become much weaker or disappear and very pronounced changes of Et 
can be expected. Having more reliable information on gl^/ is thus absolutely crucial for 
more definite conclusion on the importance of the short-range 3A^ forces considered in this 
paper. 

In Table rvnTl we show the expectation values calculated with the wavefunction obtained 
from the Hamiltonian with all 3A^ potentials in PV coupling. The expectation values of 
the TT — TT potential and of the full 3A^ potential do not vary much, but the individual 
forces show a rather strong dependence on the NN potential used. Our results (Tables 
IVII and IVIII|) for the total tt — cr and -k — uo contributions differ from those of Ref . Q| , 
where (iy(7r - a)) = 1.003 MeV and {W{n - u)) = -0.770 MeV were obtained with a wave 
function calculated with the Paris NN potential only. This difference is a consequence of the 
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AT ■not.pnfial 


Reid 


Pciris 


IVi iTriP0"PTl Q.*^ 
J. > IJ J-llCCjdJ. ty (J 


Ron n R 




—0.096 


—0.113 


—0.124 


—0.130 


7T7T(b) 


-1.822 


-1.882 


-2.086 


-1.837 


7r7r(d) 


-0.364 


-0.279 


-0.319 


-0.256 


Totcll TTTT 


-2.282 


-2.274 


-2.528 


-2.223 


/I ty I xv J I, J 






0Q7 


113 


Trr>(A+) 


0.156 


0.148 


0.164 


0.173 


/I fy\^ ) 


0.101 


0.142 


0.179 


0.174 




—0.026 


—0.036 


—0.036 


—0.037 


Total TT n 


0.287 


0.317 


0.403 


0.424 


7ro-(Z) 


0.959 


1.359 


1.632 


2.108 


TTa{N*) 


—0.497 


—0.478 


—0.581 


—0.381 


Total TTlT 


0.463 


0.881 


1.051 


1.727 




-0.902 


-1.422 


-1.722 


-2.377 




0.456 


0.453 


0.561 


0.371 


Total TTuj 


-0.446 


-0.968 


-1.162 


-2.005 


Total na + vru; 


0.016 


-0.087 


-0.110 


-0.278 


Total 3A^ potentials 


-1.979 


-2.045 


-2.236 


-2.077 



TABLE VIII: Expectation values (in MeV) of the components of the 3A^ force in Hamiltonians 
with different A^A^ potentials, calculated with eigenfunctions of the fuh Hamiltonian containing all 
listed 3A^ potentials. 



different coupling constants and cut-off parameters used in Ref. IJ]. The results collected 
in Tables IVIIII and IVll show again that also the short-range 3A^ forces are non-perturbative: 
even for the small sum of all vr — a and tt — uj contributions we get, for example with the 
Paris potential, -0.048 MeV from Table HVl and -0.087 MeV from Table ED 

Finally, in Fig.|2lwe represent the dependence of the binding energy for two different NN 
models as a function of the tcNN cut-off parameter. The short-range forces do make Et less 
cut-off dependent, compared to the case when only the vr — vr potential is included. But the 
overall effect of the n — a and n — uj potentials is rather small. We emphasize again that 
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this sum depends crucially on the values of the coupling constants used in our calculations, 
in particular on the poorly determined ratio g*/g^- 



IV. DETERMINATION OF THE LOW ENERGY CONSTANTS 
A. Effective couplings 

To connect to the 3A^ force following from ChPT, let us now consider "the pointlike 
limit", i.e., shrinking the propagators of heavy mesons to a point by taking m^AB(q3) 1 
[B = p,uj,cr), and deduce the effective contact vertices. In more detail, we separate from 
the potential two diagrams in which the pion couples to the second nucleon and re-write the 
corresponding potentials as the product 

^"(vr) {V,% + V^,) with V.^in) = ■ q2 A.(q2) , (32) 

where the factor ^'^(71) includes an overall factor i, the nNN vertex, and the pion propagator. 
Hence V^s + V^i is the vertex function of the effective point Lagrangian 

£Pomt _ (jv^FiA^) (N^TsN) Ti" . (33) 
The effective leading order Lagrangian is taken in the form 

= /:«(«) + £(^)(/5), (34) 

£(1) 

-a3 e"''^ {N^crr^N) x {N^ar^'N) ■ Vtt" , (35) 

+/?3 6"'^= vr'^ e^,, {N^a.r'N) Vk{N^a,T'N) . (36) 

It generalizes the Lagrangian considered in Ref. [l^ by including (in C^^\(3)) vertices with 

derivatives of nucleon fields.^ 

For the local part of the tt — p potential with a timelike p-exchange we thus obtain 

2 

£point^^ - pT, PV+) = , vr" (N^aN) V{N^t^N) , (37) 



^ In this paper we use the leading order ttNN couphng with the plus sign {Lt^nn — +^{N^crT°-N) ■ Vtt"), 
in agreement with Hiiber et al. 0,^3' t)ut opposite to Epelbaum et al. Therefore, the Lagrangian 

has the opposite sign compared to Eq.(2.3) of Ref. (loj. 
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Paris potential with three-nucleon force 




4 5 6 7 

A \m 1 

7[NN Jl-" 



Nijmegen 93 potential with three-nucleon force 




4 5 6 7 

A Tm 1 

7[NN Jl-" 

FIG. 2: Dependence of the triton binding energy on the TrATiV cut-off parameter, calculated for 
Hamiltonians with the Paris and the Nijmegen 93 NN potentials and various contributions to the 
3N force, which arc successively added in the following way: n — n exchange (dashed line), plus 
TT — p (dotted line) , plus ir — a (dashed-dotted line) , plus tt — u (solid line) . The actual calculations 
are indicated by the full circles, the lines are drawn to guide the eye. 
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/:p°int(^ - pS, KR) = - ^^^^ , eabc e^j t^^^N^ a^r^N) Vk{N^a,T'N) . (39) 



and comparing to ()36|) we extract 

/?2(7r-pT,PV+) = -^, (38) 

where for PS coupling there would be an additional factor 1 + Hp. For the local part of the 
7T — p KR term ()19j) one gets 

g gl{l + K p) 

For the Born ti — a exchange the point low-energy Lagrangians corresponding to the limit 
of and (El are 

2 

£point(^ - a, PV+) = ^f2_ (ATtrV N) ■ V{N^N) , (41) 

A(7r-a,PV+) = -/f^, (42) 

2 

£point^^ - a, PS-) = +77^1^ {N^N) [N^ar" N) ■ {Vn") , (43) 



/?3(7r - p5, KR) = , . (40) 



2w?m1 

ggl 



a,{n-a,PS-) = . (44) 



For TT — u exchange we obtain 



2 

^point/^ - tu, PS-) = +-^^ 71" (Art o-r'^Ar) . V(N^N) , (45) 

2m^m^ 



A(7r-c.,PS-) = (46) 

The corresponding effective interaction for the PV+ term includes an additional factor 
1 + ^ 1, which we will neglect in the following. For the contributions with an intermediate 
A^* (note that the A^* propagator is already used in static point approximation) one derives 

£P™'^V(7r- (a + cu)),N*) = ^(M^_^)(ArtAr)(Arto-r'^7V)- V7r\ (47) 



m{m* — m) ^ 

ai a + u;,N = — -( — -) . (48) 

m[m* — m) ^ ' 



It is interesting that for tt — ex exchange the result with PS tcNN coupling is closer to the 



form suggested in Refs. 



lOl isj l- Note also that the difference between the Lagrangians (PT|) 
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and ()43|) (we label it "cont" since it corresponds to the chiral contact term connecting PV 
and PS couplings) is equivalent to 

2 

£P°mt/ _ _ 99. .^tjv) [(VATt) . (jT^TT^N + ATtrW ■ (ViV)l , (49) 

This interaction is exactly the same as the second term in Q (with c — > —gg'^/{2mml.)). 
Note that the Lagrangian Q follows from the re-definition of the pion field (0), while the 
so-called chiral rotation (which transforms in lowest order from ttNN PV to PS coupling) 
is a redefinition of the nucleon field. It is argued in that the interactions in Q are 
unnatural from the point of view of ChPT and should be discarded. The appearance of such 
a term among the effective contact interactions raises the suspicion that the dynamical model 
(conforming with chiral symmetry) with a a meson and PV ttNN coupling contains also 
the mechanism corresponding to the first term in (j2)) and that, in analogy to the discussion 
in [l^, both should be discarded. Let us also point out that the most natural way of 
introducing the a meson into a chirally symmetric Lagrangian is the linear a model with 
PS ttNN coupling (which, however, fixes the coupling constant g^^). 



B. Numerical results and comparison to ChPT 

In Table IXl we list the numerical values for the dimensionless low-energy constants 

a, = AflA^ai, (50) 
A = 4^3A^A, (51) 

where = 700 MeV as in Ref. The constants following from the considered 

heavy meson exchanges appear to be of natural size. Notice that - with the exception of 
the potentials with A^* and the vr — a Born term (the latter only with PS ttNN coupling) - 
all effective vertices are of the form (jHBj) . i.e., they involve derivatives of nucleon fields. In 
particular, we do not get any contributions to the effective couplings a2 and a^. 

There is a pronounced cancellation between u and a exchanges in the low energy limit, 
both for the Born (with PV coupling) and Roper contributions. Although the u and a 
terms would always have opposite signs, their sum is rather sensitive to the poorly deter- 
mined values of the coupling constants. In particular, the almost perfect cancellation of the 
respective Roper terms appears to be accidental. The timelike part of the p exchange is 
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rather small for PV ttNN coupling, but for PS coupling the large factor 1 + Kp increases the 
corresponding (32 to the order of all the other contributions. These results are completely 
in line with our results for the corresponding potentials included in their full form into the 
Faddeev equations, i.e., the low energy limit changes of course the size of the individual 
contributions, but not their signs and relative magnitude. 

The chiral 3A^ forces and their effect on the 3A^ observables were studied in Ref. M and 

n 

more recently in Ref. [10|| . Unfortunately, the point limit of the short-range forces considered 
above cannot be compared directly to their results, because neither of these papers includes 
terms with derivatives of nucleon fields in the Lagrangian. In Ref. |]J^ it is claimed that 
these terms can be reduced to the form of (|35|) through integration by parts. To transform 
our operators into a form consistent with the effective Lagrangian of Ref. , we make the 
following momentum replacement: qs —>■ — q2. It corresponds to integrating the effective 
Lagrangians with nucleon derivatives by parts (this is equivalent to the use of momentum 
conservation at corresponding vertices, qa = — q2 — qi) and keeping only the resultant terms 
with a pion field derivative (momentum q2), while omitting at the same time still another 
term in the Lagrangian with nucleon field derivatives (leading to the momentum qi). 

However, our numerical calculations do not justify neglecting the effective Lagrangians 
with nucleon field derivatives. This is seen clearly from Table HXl where we compare the effect 
of momentum replacement in matrix elements (with the fully correlated wavefunction) of the 
original potentials, their point limits. Since these numbers are quite different, we checked 
in addition whether the disagreement is due to high momentum components of the wave 
function by multiplying the operators by the overall exponential form factor of Ref. jlfll |. 
However, we do not see any improvement with respect to the agreement between (W) and 

(W^repl). 

Nevertheless, since no other numerical estimates of LECs are currently available, we use 
the momentum replacement prescription to see how the LECs derived from the potentials 
of this paper compare with those deduced in Refs. 0, Q| . For the LECs the momentum 
replacement translates into f3i ^ a^. 

Hiiber et al. [9J introduce two LECs Ci and C2, related to our ai and as follows 

Ci = ^«i = Y^«, (52) 
C2 = ^a3 = -a3, (53) 
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Paris 




Nijmej 


;en 93 


3N potential W 


{W) 


(W,epl) 


{W) 


(W^repl) 


TT-p (KR) 


0.063 


-0.220 


0.097 


-0.231 


■K — p (KR) point limit 


0.045 


-0.223 


0.085 


-0.250 


■K — p (KR) point limit x exp. ff. 


0.053 


0.064 


0.094 


0.154 


^ - CT (Z) 


1.359 


-1.329 


1.632 


-1.617 


vr — 0" (Z) point limit 


0.873 


-0.531 


0.921 


-0.693 


vr — 0" (Z) point limit x exp. ff. 


2.034 


-0.015 


2.409 


0.186 


^ - p (T) 


-0.036 


-0.012 


-0.036 


-0.016 


TT — p (T) point limit 


-0.050 


-0.016 


-0.049 


-0.020 


TT — p (T) point limit x exp. ff. 


-0.027 


0.004 


-0.024 


0.011 



TABLE IX: Validity of the momentum replacement prescription. Matrix elements (in MeV) of 3N 
forces related by the momentum replacement prescription described in the text are compared. For 
each of the three considered cases also the point limit, as well as the point limit multiplied with 
an additional exponential form factor are shown. The matrix elements are calculated with a wave 
function corresponding the Paris NN potential together with the 3A^ forces of Table IVTm 

with A = 1 GeV. They fit these constants to Ay from nd scattering, arguing that ci is fixed to 
about ci ~ —3 (from Ay at 3MeV) and C2 ~ 0.5 . . . 1.0 (from Ay at lOMeV). From the long- 
short range potentials considered in this paper, only the KR it — p potential contributes in 
the low-energy limit to the constant C2, yielding C2 = —1.55 (see Table IX|) which differs form 
the value cited above both in sign and magnitude. However, it should be stressed once more 
that the momentum replacement was used in deriving C2{KR) to convert the corresponding 
point-like vertex from the form (j36|) to the form ()35|1 . The total Ci from Table 1x1 is Ci ~ 1, 
again quite different from that of Ref. j^. Here it may be worth mentioning that in Ref. 
the second term in the Lagrangian proportional to 02 is not considered. 

From this comparison of LECs it is at the moment hard to see whether the potentials 
used in this paper could provide a reasonable description of Ay (and if there is therefore a 
reason to believe that we have not missed some other important short-range effects). The 
comparison of LECs is seriously hampered by the fact that the ChPT derivation of Refs. 



30 





LECs 




p - T(PV) 


(32= 0.29 




p-T(PS) 


0'?= 2.17 




_ S{K¥(\ 

jj \ »- / 


3i— -2 17 


/-'nfrprill — -1 pSFi 


uj — Both 


j3i= 8.32 


ci(repl)= 2.97 


a - Born{PV) 


/3i= -6.29 


ci(repl)= -2.25 


a - Born{PS) 


ai=-6.29 


ci = -2.25 


uj-N* 


ai = 2.99 


ci = 1.07 


a-N* 


di =-2.26 


ci =-0.81 



TABLE X: Dimensionless low-energy constants defined by Eqs. (|50I51|) . derived from n — p, tt — ui 
and TT — a exchanges using the masses and coupling constants of Table (|III|) . The constants Cj(repl) 
were obtained with the help of the momentum replacement, for a — Born term this replacement is 
needed only if PV irNN coupling is used 

does not consider effective operators to which the meson-exchange potentials reduce 
to in the low-energy limit. 

In Ref. jl^ it is further argued that the Lagrangian (p!Kjl reduces effectively to just one 
independent term when used with 3A^ wave functions which are antisymmetric in the spin 
and isospin subspace of the two nucleons coupled to the contact vertex. Hence only one 
effective coupling constant (called D) for the ttNNNN vertex is considered in their numerical 
analysis. Assuming this, the matrix element of the vertex function corresponding to the 
Lagrangian is reduced to the matrix element of only one particular combination of 
spin-isospin operators. Indeed, for 

<3 = ["1 (^r^-l + ^SC^S) + "2 (TiVs + TgVi) + 2^3 (ti X T^yicT^ X CT3)] , (54) 

it holds that 

= »l-»2^+^»3 ^^^^ _ _ _ ^ ^ ^^^J ^ ^55^ 

where Ai^ ~ ~ -^13-^13) spin-isospin antisymmetrization operator for particles 1 

and 3. Therefore, if only such components of 3A^ wave functions are important, the effect 
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of ()35p can be represented by a single term. However, we would like to point out that the 
3A^ wave function contains also components which are symmetric in spin-isospin space (and 
antisymmetric in their spatial part). For such components one gets instead of (f^ 

_^a^-a^-Aas ^^^^ _ _ ^ ^ ^ ^^^^ ^ ^^^^ 

where 5*13 = |(1 + PfgP^^) is the spin-isospin symmetrization operator for particles 1 and 3. 
If these components are not omitted, the matrix elements of the potentials derived from the 
ttNNNN interaction ()35|1 still contain three different combinations of the constants ai,a2, 
and as. The components of the triton wave function with odd orbital angular momentum 
between nucleons 1 and 3 contribute about 5 per cent to the total norm. Even though the 
weight of these states in the norm is not large, their effect can become significant in matrix 
elements of operators with the proper tensor structure. Moreover, NN P-waves are known 
to very important for the description of Ay in Nd scattering. 

If we neglect for the moment the symmetric spin-isospin wave function components and 
interactions with nucleon field derivatives, then the effect of the terms is reduced to the 
calculation of the matrix element of In this approximation, the dimensionless constant 
cd oi Ref. [l^ is given in terms of a, as 

cd = -ai + a2- 4a3 . (57) 

From Table [X] we get cd ~ 6 ... 8, depending on the type of the ttNN coupling. This 
numerical value of co is larger than cu ~ 1.8 . . .3.6 obtained in Ref. In fact, the KR 
contribution dominates our value of cd, due to the factor 4 in the last term of ()57|) . As in the 
comparisons above to Ref. we had to compare quite different effective Lagrangians and 
the approximations made in establishing the correspondence between the effective constants 
appear to be too crude to allow drawing clear conclusions. 

Finally we note that the efforts of determining the LECs in may have resulted in a 
better description of Ay with different constants, had the additional freedom in the param- 
eters of the TM n — n force due to the experimental uncertainties discussed in Sec. IIII AI 
been taken into account. 
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V. CONCLUSIONS 



The conclusions of this paper can be summarized in three points: 

1. When considering the effect of the TM n — n exchange 3N force on the triton binding 
energy, only the variation of the cut-off parameter, which is often adjusted to reproduce the 
experimental value, is studied in the literature. 

However, the tt — tt force contains the vrN scattering amplitude as a building block, which 
at low pion momenta has a model independent form parameterized by three constants which 
have to be extracted from experiment. In contrast to the cut-off parameter, each of these 
three constants multiply different spin-isospin operators, and therefore act differently on the 
various channels of the wavefunction. We calculate the propagation of the experimental 
errors of the nN data, which are used to extract the off-shell TM nN scattering amplitude, 
into the 3N force parameters. As a consequence, the triton binding energy calculated with 
the TM TT — TT potential has an uncertainty of about ±0.4 MeV, which is almost entirely due 
to the experimental errors in the nucleon sigma term. 

2. The long-short range n — a and tt — a; exchange 3A^ forces individually have large 
effects on the triton binding, but two kinds of cancellations determine their net effect. We 
find contributions of opposite sign and comparable or very close in magnitude from Z-graph 
and A^* excitation, respectively, in the ir — a and n — uj 3N forces. Also Z-graph and A^* 
potentials cancel in part for each of these long-short range forces separately. 

The extent of the cancellation is controlled decisively by the ratio between the 
cr7V*(1440)A" and a;A"*(1440)A" couphng constants. If this ratio is taken to be the same 
as for the aNN and uNN couphng constants, as suggested by the naive constituent quark 
model, the attraction of the uj exchange, originated by the Born terms of relativistic order, 
dominates. 

The change in binding energy caused by these 3A^ potentials cannot be calculated pertur- 
batively, since their inclusion in the Hamiltonian changes the wavefunction significantly. For 
instance, we find that, although the net effect of tt — o" and tt — cu exchange 3A^ forces on the 
binding energy is small, the resulting wavefunction yields rather different matrix elements 
of individual interactions from the ones calculated without those terms in the Hamiltonian. 
Therefore, it would be very interesting to study these variations in the wavefunction, for 
instance in electron scattering on the 3A^ bound state. Furthermore, these potentials may 
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strongly influence low-energy Nd scattering observables. 

In the case of PS ttNN coupling, both Z-diagrams are attractive, resulting in a strong 
overbinding of the triton. In the future, these Z-diagrams should be complemented by 
additional potentials derived from chiral contact ttNNcx and ttNNuj vertices. 

3. We extracted the LECs from the low-energy limit of the long-short range meson 
exchange potentials and compared with the ChPT results of Refs. jo], |lO|. The constants 
we obtain are of natural size. Unfortunately, a direct comparison with the ChPT results is 
not possible since we arrive at contact vertices of different structure. To translate them into 
the ChPT form certain approximations are necessary, the validity of which we found to be 
unfounded in our numerical evaluation. It does not come as a surprise that our LECs do not 
agree with those of Refs. Efforts should be joined from the two sides (meson-exchange 

description and ChPT methods) to investigate this question further. 



Acknowledgments 

J.A. was supported by the grants GA CR 202/00/1669, GA CR 202/03/0210 and by the 
Lisbon project. He would like to thank his Portuguese colleagues for the warm hospitability 
during his stay in Lisbon. M.T.P. and A.S. were supported by FCT and FEDER under 
grants CERN/FIS/43709/2001 and POCTI/FNU/40834/2001. 



APPENDIX A: CONTRIBUTIONS OF NUCLEON BORN DIAGRAMS 



We will give here some details of the derivation of nucleon Born diagram contributions to 
the 3A^ force from the exchange of one pion and one heavy meson (B), since we experienced 
some misunderstandings and many questions when discussing their origin. Our treatment 
is a condensed version of the technique developed in a number of papers on e.m. meson 



exchange currents 
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36, 



We pay special attention to the dependence of our results on 
the type of the nNN vertex. 

The generic Feynman amplitude W(l) corresponding to Fig. reads: 



>V(1) = J^''P'' + 2^3, 



(Al) 



where we gather all factor connected with the "active" first nucleon into the amplitude 
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and the rest into the amphtude T> 

= - 7i(p'2)r (tt, g2)n(p2) u{P3)^iB, qsHps) ^^2) ^sigs) , (A2) 
^'^ = u{p[)[T'^{7r,-q2) S{P')T{B,-q,) + T{B,-q,) S{P)T'^i7T,-q2)]u{p,) , (A3) 

where S{p) = {i'y^p^ + m)~^ and Asiq) = {rn\ + q^)~^ (the Pauh metric is used), P = 
Pi + = Pi — 12, P' = Pi — = Pi + <i2 with qi = p[ —pi, and we suppress (unless we specify 
B) additional isospin and Lorentz indices, that would appear for non-scalar and charged 
meson B. The pion {ttNN) vertex functions (corresponding to iC) are 

r(7r-PS) = -^?75r^ (A4) 
r'^(7r-PV,g) = ^^g^7^75r% (A5) 

where q = p' — p is the momentum of the pion entering the vertex (the signs of our vertex 
functions differ from those of Ref . js^l , we adopted the current convention to agree with that 
usually used in ChPT calculations). The vertex functions r{B,q) for the heavy mesons B 
will be specified later. Note, that from the identities 

^(p;)r'^(7r-PV, -q^) S{P') = ^(p;)[r«(7r-PS) S{P') + £^Tor^] , (A6) 

S{P) r"(7r-PV, -g2)«(Pi) = [S{P) T^{n-FS) + £-^,T^]u{pi) , (A7) 

it follows 

J-'py = JFpg + ^cont 5 (A8) 

^cont = ^^(P'l) {75^r, r(i?, -g3)}^«(pi) . (A9) 

The "contact" amplitude arises from the usual chiral contact interaction, which appears in 
the chiral rotation from the PV to PV ttNN coupling. 

To get the quantum- mechanical 3A^ potential from the Feynman amplitude ()A1|) . it is 
necessary to subtract the part which is in the quantum mechanical description already 
included in the T matrix in the iterations of the A^A^ potential. We split the nucleon 
propagator into its positive and negative energy parts S{p) = S^{p) + S^{p): 

^ ~ HfiPfj. ^ + -^74 — «7 ■ p rn — £'74 — ^7 ■ p u{p)u{p) v{—p)v{—p) 



Sip) 



p^ + m^ 2E{E-po) 2E{E + po) E - po E + po ' 

(AlO) 
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(where E"^ = p^ + m^) and define the corresponding amplitudes W^(l) and J-'^. The spinors 
u and V are given by: 

For the calculations of this paper it is sufficient to keep only the leading order in p/m and 
replace E ^ m. 

The true "pair" (or "Z-diagram" ) contributions to the 3A^ potential are defined by the 
straightforward non-relativistic reduction (i.e. the decomposition of the spinor matrix el- 
ements in powers of p/m keeping only leading order terms) of the VV^(l). In the order 
considered, this contribution is non-zero only for the PS ttNN coupling for any meson ex- 
change. For this coupling (and again, at the given order in p/m) VV^(l) with positive 
energy nucleon in the intermediate state corresponds to the iteration of the lowest order 
non-relativistic OPEP potential (which is independent of the type of ttNN coupling). On 
the other hand, for the PV ttNN coupling, the pair diagrams VV^(l) do not contribute. But 
the positive energy part VV^(l) does not exactly equal to the iteration of the non-relativistic 
OPEP, since in the Feynman amplitude the vertex function is off-mass-shell (P^ 7^ 7^ P'^) 
while in the iteration of the OPEP this potential is off-energy-shell. One can re-arrange the 
energies in ^20 entering the PV vertex function in ()A1IA3|) identically 

g^o = P^- E{p[) = {E{P') - E{p[)) - {E{P') - P^) , (A12) 
g2o = E{p,) - Po = (P(pi) - P(P)) + (P(P) - Po) . (A13) 

The first energy differences on the r.h.s. put the PV vertex on its mass shell. The corre- 
sponding part of W'''(l) is then again identified with the iteration of the non-relativistic 
OPEP and it is identical to the full W"'"(l) for the PS coupling. But the second terms on g2o 
cancel the denominator of S^{P') or S^{P) and give rise to a contribution to the quantum 
mechanical potential Vr^(l). 

To sum it up, the relevant contributions to the 3N potential are 

W{1) = E" 0^ + 2 ^3, (A14) 
^ £1 ^2 (^^2 ■ q2) w(P3)r(5, g3)w(P3) Kiq2) ^BiQs) , (A15) 
F^k ^ £[^(p'i)75r>(-P')^^(-P')r(P,-g3)«(pi) + 
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u{p[)T{B, -qsH-P) t;(-P)75rXpi)] , (A16) 

i^pV ^ -^[^(p'i)7475r>(P')^(P0r(5,-g3)«(pi)- 

u{p[)T{B, -qs)u{P) ^i(P)7475rrn(pi)] , (A17) 

^cont ^ ^^(P'i){^5rf, r(5,-g3)}^«(pi). (A18) 

where ~ stands for the non-relativistic reduction. Since only F*^'" survives for PS couphng 
and only contributes for the PV one (and since the iterations of non-relativistic OPEP 
are identical in both cases), one gets from ()A9jl 

FpV = FpV + F:,,, , (A19) 

and therefore 

W^P+v = + W,ont , (A20) 

for all heavy meson exchanges considered. 

It remains to list the vertex functions for the a, u and p mesons and the results of the 
non-relativistic reduction. For the scalar isoscalar a meson one gets 

T{a) = tg^, (A21) 

D'^ia) = T^{cT2 ■ q2) A.(g2)A<,(g3) , (A22) 

i^pVW = -^^^i"('^i-q3), (A23) 

Fphi^) = (A24) 

Frontier) = +^^n"(^i-qi), (A25) 

where ()A19j) can be verified with the help of qi + q2 + Qs = 0. The corresponding potentials 
are just the products of D and F factors, they are listed in the main body of the paper. 

For the isoscalar vector meson u only the timelike part /i = 4 of the vertices contribute 
up to the order considered 

r/.(^,g) = gu,[l^, - ^^M^^i'] ' (A26) 
Dliuj) = 1^ T^{cT2 ■ q2) A.(g2)A^(g3) , (A27) 
FpU^) = -^^f('^i-q3), (A28) 

^Ps:4(^) = -^-^^ r^i^i ■ qs) , (A29) 
= +^^i"(^i-q3). (A30) 
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Finally, for the isovector vector meson p the vertex function reads 

ri(p, q) = 9p[l, - ^^^^A ■ (A31) 
For the timelike component of the p field [p = 4) one gets 

Dtip) = II T^r'.ia, ■ q,) A,{q,)A,{q,) , (A32) 

Fp'^AiP) - [^-^^""^ ■ - Wt{<ri ■ Qi)] , (A33) 

F^Up) = -^^"1^ ■ ^3) , (A34) 

^cont,4(p) = +^ [<^a6«p(0-l • Qa) + ^eaOcTf (tTi • Qi)] , (A35) 

where Qi — p- + pi. For the spacelike component 

D'^^p) = T!^Tl{cr^ ■ qs) [Qs + (1 + ^p)i<Tz x qs] A,(g2)A,(g3) , (A36) 

Y%^{p) = 0, (A37) 

Ytf{p) = -F:L(p) = e^^cTicT, . (A38) 



In the main body of the paper we denote the contributions due to the forth component of 
vector fields by subscript T and those from the exchanges with // = 1, 2, 3 by the subscript 
S. 



APPENDIX B: cr-MESON EXCHANGE IN 3N POTENTIALS 

The Bonn potentials use different parameters (masses, coupling constants and cut-off 
parameters) of their "cr" -meson in NN channels with isospin 7 = and I — 1. Let us define 



9l(I) F^(I)NNi^'^{I)NN , q^) 

47r ml^j^ + q2 



v{a{I), q^) = ^ ^^^^^^^7 ^..--^ ^ (g^) 



where ga-{i), Af^(^j)j^r^,'m'^^^ are the Bonn (t( 7) -exchange parameters in the respective NN 
isospin channel. The dominant central part of the Bonn cr-exchange potential is then given 

by 

vcici') = v(a(0), q^) Po + v{a(l), q^) P, , (B2) 
Po = \{l + T,-T2), Pi = i(3-Ti-r2), (B3) 
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FIG. 3: The fit (solid line) to the central part of the isoscalar Bonn potential (short-dashed line; 
practically coincides with the solid line on the graph), the first term of Eq. (|B4|) . For comparison, 
also v{a{0),q^) (dash-dotted line), v{a{l),q^) (dashed line) and their isovector combination (the 
second term in Eq. HB4|) : dotted line) are shown. 

which can be re- written in a form similar to the central potentials originated from exchanges 
of two scalar-isoscalar and two scalar-isovector particles: 



This does not mean that the Bonn a-exchange is equivalent to the true exchange of four 
scalar particles, since: 1) the last term has the sign opposite to the sign of a real scalar 
exchange; 2) all four of these exchange have to act at once, i.e., they cannot be separated 
by an exchange of another meson (in iterations of the Lippmann-Schwinger equation). 

To have a simple prescription for the a-exchange in our 3A^ potentials, we fitted the first 
term in Eq. ()B4j) by the function v{a, q^) dependent on parameters of single true isoscalar 
cr-exchange g^, AaNN,ml. The result of the fit to the Bonn B potential gives the values 
listed in Table IIIH and the quality of the fit is shown on Fig. El The fitted parameters 
are between the Bonn values for NN / = 0, 1 channels. The second ( "isovector-exchange" ) 




1 r 



'v{a{0), q2) + 3v{a{l), q^)] + n ■ T2 ^ [t;(a(0), q^) - v{a{l), q^)] . (B4) 
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component of Eq. ()B4|) is much smaller in absolute value and even negative for q < 10/m~^ 
(see Fig. E)), which precludes an approximation by the exchange of scalar particle(s). If we 
nevertheless include it into the 3A^ potential, it gives a considerably smaller contribution to 
the triton binding energy. Therefore, we neglected it in the calculations of this paper. 
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